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ABSTRACT 

When a gravitationally lensed source crosses a caustic, a pair of images is created or destroyed. We calculate 
the mean number of such pairs of micro-images (n) for a given macro-image of a gravitationally lensed point 
source, due to microlensing by the stars of the lensing galaxy. This quantity was calculated by Wambsganss, 
Witt & Schneider (1992) for the case of zero external shear, 7 = 0, at the location of the macro-image. Since in 
realistic lens models a non-zero shear is expected to be induced by the lensing galaxy, we extend this calculation 
to a general value of 7. We find a complex behavior of (n) as a function of 7 and the normalized surface mass 
density in stars k*. Specifically, we find that at high magnifications, where the average total magnification of the 
macro-image is = |(1-k*) 2 -7 2 | _i ^> 1, (n) becomes correspondingly large, and is proportional to The 
ratio (n) /(/i) is largest near the line 7 = 1 - k» where the magnification becomes infinite, and its maximal 
value is 0.306. We compare our semi-analytic results for («) to the results of numerical simulations and find 
good agreement. We find that the probability distribution for the number of extra micro-image pairs is reasonably 
described by a Poisson distribution with a mean value of (n), and that the width of the macro-image magnification 
distribution tends to be largest for (n) ~ 1. 

Subject headings: cosmology: gravitational lensing, dark matter — quasars: general 



1. INTRODUCTION 

Gravitational microlensing by the stars of a lensing galaxy 
can have a large effect on the magnification of lensed sources 
I (Chang & Refsdal 1979; Young 1981; Paczyriski 1986). As the 
macro-images of multiply imaged sources are typically located 
in relatively dense star fields of the lensing galaxy, microlensing 
is quite common in such systems. The typical angular separa- 
tion between the micro-images of a cosmological source due to 
a stellar mass micro-lens is of order 1 micro-arc-second, which 
, is far too small to be resolved. For the time being, the only 
observable manifestation of microlensing is to change the mag- 
nification of the macro-image relative to the average magnifica- 
tion that is predicted for the smoothed out surface mass density 
profile of the galaxy. 

The first observational evidence for quasar microlensing 
was found by Irwin et al. (1989) in the quadruple system 
Q2237+0305, which subsequently has been monitored by many 
groups (Corrigan et al. 1991; Burud et al. 1998; Lewis et al. 
1998, Wozniak et al. 2000a,b). In particular the latest results 
show that all four quasar images vary dramatically, going up 
and down by more than one magnitude on timescales of less 
than a year. The fact that individual (caustic-crossing) events 
can be clearly distinguished allows to put upper limits on the 
source size (Wambsganss, Paczyhski & Schneider 1990; Yone- 
hara 1999, 2001; Wyithe et al. 2000). 

In the double quasar Q0957+561, originally there was an al- 
most linear change detected in the (time-shifted) brightness ra- 
tio between the two images (Aotab ~ 0.25 mag over 5 years), 
which was interpreted as microlensing by solar type stars. But 
since about 1991, this ratio stayed more or less "constant" 
within about 0.05 mag, so not much microlensing was going 
on in this system recently (Schild 1996; Pelt et al. 1998). How- 



ever, even the "lack of microlensing" in this system can be used 
to put limits on compact dark matter in the halo of the lensing 
galaxy (Wambsganss et al. 2000). 

A number of other multiple quasar systems are being mon- 
itored more or less regularly, with some showing indications 
of microlensing, e.g. H1413+117 (0stensen et al. 1997), 
B0218+357 (Jackson et al. 2000), or HE1 104-1805 (Gil- 
Merino et al. 2002, Schechter et al. 2002). For a recent review 
on quasar microlensing see Wambsganss (2001). 

Microlensing has recently been suggested as the source of 
short time scale low-level variability in the time delayed light 
curves of multiple images of quasars (Wyithe & Loeb 2002), 
and might help study the properties of broad-line clouds in 
quasars. A better understanding of the microlensing by stars in 
galaxies would provide a better handle on the probability dis- 
tribution for microlensing of cosmological sources (Wyithe & 
Turner 2002) that is relevant for Gamma-Ray Bursts (GRBs), 
especially in light of a possible microlensing event that was ob- 
served in the optical and near IR light curve of the afterglow of 
GRB 00030 1C (Garnavich, Loeb & Stanek 2000; Gaudi, Gra- 
not & Loeb 2001; Koopmans & Wambsganss 2001) 

Another example where microlensing may play an impor- 
tant role is in explaining the flux ratio anomalies observed in 
close pairs of images of quadruply lensed quasars (Schechter & 
Wambsganss 2002). Such systems are usually modeled using a 
simple smooth surface mass density profile for the galaxy, pos- 
sibly with the addition of an external shear. While these mod- 
els successfully reproduce the observed locations of the macro- 
images, the flux ratios they predict are quite often in poor agree- 
ment with observations. Specifically, the theoretical flux ratio 
for a close pair of highly magnified macro-images is 1:1, while 
observations show a difference of up to one magnitude. For ex- 
ample, MG0414+0534 has an observed flux ratio of 2:1 in the 
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optical (Hewitt et al. 1992; Schechter & Moore 1993), while 
in the radio the flux ratio is 1:1 (Trotter et al. 2000). There 
are also alternative explanations for these flux ratio anomalies, 
such as intervening dust (Lawrence et al. 1995) or milli-lensing 
by galactic sub-structure (Mao & Schneider 1998; Metcalf & 
Madau 2001; Dalai & Kochanek 2002; Chiba 2002). The 
study of microlensing can help distinguish between these al- 
ternative explanations, and may be useful in constraining the 
sub-structure of galaxies. 

In many cases the study of microlensing cannot be done an- 
alytically, and much of the work is done using numerical sim- 
ulations. The magnification distributions of the macro-image 
(MDMs) are particularly important for understanding the ob- 
served properties of lensed systems. The latter have been cal- 
culated for a wide range of parameters using numerical simu- 
lations (Wambsganss 1992; Rauch et al. 1992; Lewis & Irwin 
1995, 1996; Schechter & Wambsganss 2002). While simula- 
tions are applicable to a wide range of problems, and are in 
many cases the only available technique, they are usually time 
consuming and do not always provide a qualitative understand- 
ing of the results. Specifically, they do not seem to provide 
an explanation for the detailed structure that is present in the 
MDMs that they produce. 

In a recent paper, Schechter and Wambsganss (2002) ex- 
plained the flux ratio anomalies as resulting from a differ- 
ent qualitative behavior of the MDMs for macro-minima and 
macro-saddlepoints in the arrival time surface, that results in a 
larger probability for de-magnification (relative to the average 
magnification) of saddlepoints, compared to minima. This is 
in agreement with the observations of MG0414+0534 and four 
other recently discovered quadruple systems (Reimers et al. 
2002; Inada et al. 2002; Buries et al, in preparation; Schechter 
et al., in preparation) in all of which the fainter image in a pair 
corresponds to a saddlepoint, while the brighter image corre- 
sponds to a minimum. 

The structure of the MDMs seems to be tightly related to the 
probability distribution for the number of extra micro-image 
pairs (EIPs) corresponding to a given macro-image (Rauch et 
al. 1992). This may be expected since the magnification of the 
macro-image is simply the sum of the magnifications of all the 
micro-images that it is composed of. An analytic expression for 
the mean number of EIPs was derived by Wambsganss, Witt & 
Schneider (1992, hereafter WWS) for the case of zero external 
shear. In §2 we generalize this result for arbitrary values of the 
external shear, and obtain a semi-analytic expression. More de- 
tailed expressions are provided in Appendix A, along with an 
analytic result for the case of equal shear and convergence in 
stars. In §3 we compare our analytic results to numerical simu- 
lations, and find good agreement. We also show that the prob- 
ability distribution for the number of extra micro-image pairs 
is reasonably described by a Poisson distribution. The possible 
implications of our results are discussed in §4. 

2. ANALYTIC RESULTS 

In this section we calculate the mean number of positive par- 
ity micro-images (N), and the mean number of extra micro- 
image pairs (EIPs) («), that are induced by the random star field 
of the lensing galaxy near the location of a macro-image. The 
macro images are located at stationary points (i.e. minima, sad- 
dlepoints or maxima) of the time delay (Fermat) surface of the 
smoothed out surface mass density distribution of the lensing 
galaxy. The source size is assumed to be small compared to 
the Einstein radius of the stars in the galaxy and compared to 



the typical distance between caustics, so that we may assume 
a point source. In this limit the probability distribution of the 
random shear caused by the star field and the resulting distribu- 
tions of the number of micro-images and total magnification of 
the macro-image are independent of the mass spectrum of the 
stars (Schneider & Weiss 1988). The remaining parameters on 
which the microlensing characteristics in this case might still 
depend are the normalized surface mass density in stars k* and 
in a smooth component (dark matter) well as the large 

scale shear of the galaxy 7. However, Paczyhski (1986) has 
shown that such a model is always equivalent to a model with 
effective convergence and shear given by 



K e Jf = 



(1-Kc) 



and with no smooth component {n e / f = 0). The effective mag- 
nification in this model is related to the true magnification by 
[fff = ^(1 - k c ) 2 . Therefore, without loss of generality, we re- 
strict ourselves to k c = 0, while letting and 7 vary. 

Each pair of micro-images (EIP) consists of a saddlepoint 
and a minimum, and thereby includes one positive parity image. 
For a macro-minimum there is also an additional positive parity 
micro-image corresponding to the global minimum. Therefore, 
n=N- \ for a macro-minimum and n = N otherwise, where n 
is the number of EIPs, and N is the number of positive parity 
micro-images. This also carries through to the average values 
of these quantities: 
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(N)-l k* + 7<1 
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K» + 7 > l 



It is therefore sufficient to calculate one of these quantities in 
order to determine the other. The values of (N) and (n) for the 
case of zero external shear (7 = 0) were calculated by WWS. In 
this work we follow their analysis and generalize their result to 
the case of a non-zero shear 7. The average total magnification 
of the macro-image is given by 
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|(l-«*) 2 -7 2 | 



(3) 



As the magnification is solely due to area distortion, flux con- 
servation implies that a sufficiently large area Ad in the deflec- 
tor plane will be (backward) mapped onto an area A s = Ad / (/x) 
in the source plane (when projected onto the deflector plane, 
so that it extends the same solid angle from the observer). The 
mean number of positive parity micro-images (N) is equal to 
the multiplicity q by which regions of positive parity within 
Ad, when mapped onto A s , cover A s . 

The probability distribution of the random shear produced by 
the stars is given by 



p(K*,S u S 2 ) = 



27T( K l + S 2 l + 52)V2 



(4) 



(Nityananda & Ostriker 1984; Schneider, Ehlers & Falco 1992) 
where S\ and S2 are the two components of the internal shear. 
As we assume that k c = and that the stars are point masses, 
the convergence vanishes everywhere (except at the locations 
of the stars, where it is infinite, but these form a set of mea- 
sure zero for any finite area on the deflector plane) and the local 
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magnification at a given location in the star field is just due to 
the shear at that point 



M7>Si,S 2 ) = 



1 



l-(7+Si) 2 



si 



(5) 



where for convenience we have chosen Si to lie in the direction 
of the external shear 7. The area dad = AdP(K*,S\,S2)dS\dS2 
within Ad where the shear lies between (7+51,^2) and (7+ Si + 
dSi,S2 + dS2), is mapped onto the area da s = dad/ ' P(l,Si,S2)- 
Therefore 



2tt|(1 



p(n*,S\,S2)dS\dS2 



M>0 



(6) 
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l-( 7 +5i) 2 -S 



(«2+S 2 + 5 2 )3/2 



This integral may be evaluated analytically. However, the re- 
sulting expression is long and cumbersome so that we prefer not 
to write it down explicitly. Instead, we show contour plots of 
(«) and (n) / (p) in Figure 1, and provide the values of (N) / (p) 
at representative points of (ac*,7) in Table 1. In Appendix A 
we reduce the expression in equation (6) to a one dimensional 
integral in two different ways (so that (N) may be easily eval- 
uated numerically) and provide an analytic expression for the 
case where 7 = ac* . For 7 = 0, equation (6) reduces to a simple 
analytic expression (WWS): 



(AO 



= (AOU-k*) 2 = 1+2k 2 -2k» v / / 1 + K 2 . 



(7) 



For = 0, there are no extra images due to stars, and we have 
either one positive parity image for a macro-minimum (N = 
(N) = 1 for 7 < 1) or none for a macro-saddlepoint (N = (N) = 
for 7 > 1). Hence, for K,=0we have: 
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7 < 1 
7> 1 



(8) 



As can be seen from equation (6), the ratio (N) / (p) remains 
finite and varies smoothly with 7 and ac* near the lines of in- 
finite average total magnification (p) in the 7-Ac* plane (i.e. 
7 = 1 1 — k* |). Together with equation (2) this implies the same 
for the ratio (n)/(p) along the line 7 = /c* - 1, while along the 
line 7 = 1 — ac* it is continuous but its derivative is discontinuous 
in any direction that is not along this line (as can be seen in the 
lower panel of figure 1). 

Furthermore, (n)/(p) decreases on either side of the line 7 = 
1 - ac*, and therefore attains its maximal value along this line, 
((n)/(/i)) max = 0.306136 at (ac*, 7 ) = (0.37895,0.62105). For 
(p) » 1 (7 w |1-ac*|), small changes in7orAc* can cause large 
changes in (p) and (n), while the ratios («) / (p) or (N) / (p) re- 
main approximately constant. In this region, when crossing the 
line of infinite magnification 7 = 1 - ac*, from a macro-minimum 
(7 < 1 - ac*) to a macro- saddlepoint (7 > 1 - ac*), (N) decreases 
by 1 (corresponding to the macro-minimum that disappears), 
but since (TV) is infinite at this line, this amounts to a zero frac- 
tional change in (N) . 

3. COMPARISON TO NUMERICAL SIMULATIONS 

In this section we compare the analytic results of §2 to the re- 
sults of numerical simulations. Combining a ray shooting code 



(Wambsganss 1990, 1999) with a program that detects the lo- 
cation of the caustics (Witt 1993) as in WWS, we extract ad- 
ditional information on image multiplicity and magnification. 
Whenever the source crosses a caustic, a pair of images consist- 
ing of a micro-saddlepoint and a micro-minimum are created 
or destroyed. Each simulation is based on a particular realiza- 
tion of the random star field, that determines a caustic network 
in the source plane, and in principle determines the number of 
such extra image pairs (EIPs) at any point in the source plane. 
In order to calculate the average number of EIPs, (n), from the 
results of a simulation, we identify the regions with different 
number of EIPs, n, on the source plane, and calculate the frac- 
tion /„ of the source plane that they cover, which is equal to 
the probability p„ of having n EIPs at a random location on the 
source plane. This is illustrated in Figure 2. Pixels in the source 
plane that are crossed by caustics (which are colored in yellow 
in Figure 2) are attributed to the corresponding higher image 
number. There are also parts of the source plane for which it is 
very difficult to uniquely identify n (corresponding to the black 
regions in Figure 2) due to the occasionally complex caustic 
structure, combined with the finite pixel size of the simulation. 
These regions are left out when calculating /„. The "uniden- 
tified" regions correspond to <~ 1% of the source plane, and 
a large fraction of these regions probably corresponds to rela- 
tively high values of n. However, as we do not know exactly 
what distribution of multiplicity n we should assign to these re- 
gions, and lacking a better option, we simply assign to them the 
same distribution of the identified regions: 



Pn = 



fn 



(9) 



This is equivalent to leaving out the unidentified regions en- 
tirely, and normalizing the probability distribution of p n . The 
average number of pairs for the simulation is then given by 



(»)sim = y^ftj?n 



(10) 



Since we have simple analytic expressions for (n) for 7 = 
(WWS, see equation 7) and for ac* = (equation 8) we chose 
values of (ac*,7) along the line ac* = 7 for the simulations, so 
that they would serve as a good check for our analytic re- 
sults. An additional advantage of this choice is that it corre- 
sponds to the interesting case of a singular isothermal sphere 
(for ac c = 0). We performed three simulations, two for macro- 
minima (ac* = 7 = 0.333, 0.400 with (p) =3, 5) and one for 
a macro-saddlepoint (ac* = 7 = 0.666 with (p) = 3). The re- 
sults of the simulations are shown in Tables 2 and 3. For 
ac* = 7 = 0.333,0.400 and 0.666, (n) sim is 3.7%, 3.3% and 
7.7%, respectively, lower than the analytic result (n) from §2 
(denoted by (n) t h in Table 2). The values of («) s im in Table 2 
were calculated according to equation (10) which assigns the 
p n distribution of the regions with identified image multiplicity 
n in the source plane to the unidentified regions. However, as 
mentioned above, the unidentified regions are typically related 
to complicated caustic structures, and are hence more likely to 
contribute to larger values of n compared to the identified re- 
gions. If the average n of the unidentified regions is say, 5, 
this would make s im 1.1% higher, 2.6% higher and 1.1% 
lower than (n), for ac* = 7 = 0.333, 0.400 and 0.666, respec- 
tively. The latter should assume values of 3.96, 3.19 and 5.70, 
for ac* = 7 = 0.333, 0.400 and 0.666, respectively, in order for 
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(n)sim to be exactly equal to («). One should also keep in mind 
the "cosmic" variance, i.e. the fluctuations between the results 
of different simulations for the same and 7, due to different 
statistical realizations of the star field over a finite region in the 
deflector plane. We therefore conclude that there is good agree- 
ment between the results of the numerical simulations and our 
analytic results for the average number of EIPs (n) . 

The average total magnifications from the simulations, 
(/x)sim, are 3.1%, 0.16% and 3.5% lower than their theoretical 
values for k* = 7 = 0.333, 0.400 and 0.666, respectively. The 
scatter in (fi) S im can be attributed to the "cosmic" variance. The 
fact that (/i) S im is on average slightly lower than its theoreti- 
cal value arises since we consider finite regions in the deflector 
plane and in the source plane. Rays that fall very close to a star 
in our deflector field suffer very large deflection angles which 
may take them outside of our source field; these are not com- 
pensated for by rays with large deflection angles from stars out- 
side of our deflector field, that should have been deflected into 
our source field (see Katz, Balbus, Paczyhski 1986; Schneider 
& Weiss 1987). We conclude that the numerical simulations are 
in good agreement with the theory on the value of (p), as well. 

The magnification distribution of the macro-image (MDM), 
pip), where p(p)dn is the probability that the total magnifi- 
cation of the macro-image is between [i and fj, + dp, can be 
expressed as a sum over the contributions from regions with 
different numbers of EIPs n, 



P(.V) = ^Pn(p) , 



(ID 



n=0 



where p„(p)dp = (dp„/dp)dp is the probability of having n 
EIPs and a total magnification between fj, and fi + dp,, with the 
normalization 
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p{p)dp=\. (12) 



In Figure 3 we show the results for p n from our simulations 
(that are given in Table 2), along with a Poisson distribution: 



Pn = 



(13) 



where the solid line is for a = (n) and the dashed line is for 
a = (n)sim, calculated according to equation (10). A Poisson 
distribution provides a reasonable fit to the results of all our 
simulations. For large values of n there is a relatively larger de- 
viation from a Poisson distribution. This results, in part, from 
the difficulty in identifying regions with large n in the source 
plane. One should also keep in mind that the "cosmic" vari- 
ance in p n becomes larger with increasing n. There still seem 
to be some systematic deviations from a Poisson distribution, 
however they appear to be small for relatively low values of n, 
which cover most of the source plane. We therefore consider a 
Poisson distribution for p„ to be a reasonable approximation. 

Figure 4 shows plots of pp(p) and pp n {p) from our simula- 
tions. The shape of p n (p) for different n seems quite similar, 
while the average magnification, 
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{p)n =— p n {p)pdp 
Pn JO 

and the normalization, p„, are different. 



(14) 



4. DISCUSSION 

We have calculated the mean number of extra micro-image 
pairs {n) of a point source, as a function of k* and 7. The 
results are shown in Figure 1 and Table 1. One dimensional 
integrals for general values of (k*,7) and an analytic result for 
the case 7 = are presented in Appendix A. Near the lines of 
infinite magnification in the 7 - ft* plane (7 = 1 1 — k* | ), (n) di- 
verges, and is proportional to the mean macro-magnification 
(p). The ratio (n)/(p) is continuous along these lines and 
varies smoothly along the line 7 = k* — 1, while its derivative 
is discontinuous along the line 7 = 1 - in directions that are 
not along this line. This creates a "ridge" in (n) / (p) along the 
line 7=1- «;*, where it also peaks at (k*,7) = (0.379,0.621) 
with ((n)/(//))max = 0.306. The analytic results for (n) are in 
good agreement with the results of numerical simulations we 
performed for k» = 7 = 0.333, 0.400 and 0.666 (as can be seen 
in Table 2). 

We find that the probability distribution p n for the number 
of extra image pairs (EIPs) n, that is calculated from numerical 
simulations, may be reasonably described by a Poisson distri- 
bution. This result holds both for the numerical simulations 
performed in this paper (e.g. Table 3 and Figures 2 and 3), and 
for numerical simulations from previous works (Rauch et al. 
1992). Furthermore, the shape of the magnification distribution 
p n {p) of regions with a given n appears to be similar for dif- 
ferent values of n, where only the overall normalization p n and 
mean magnification (p)„ depend on n (see Figure 4). 

The mean number of EIPs (n) can serve as a rough mea- 
sure for the width of p(p), the magnification distributions of the 
macro-image (MDM). For (n) <C 1 there is little contribution to 
the MDM from regions in the source plane with n > 0, since 
these regions cover only a small fraction of the source plane. 
For («) ^> 1 the Poisson distribution p„, for n, approaches a 
Gaussian distribution with a mean value of («) and a standard 
deviation of a = («) 1//2 , so that only ~ (n) 1 / 2 different values 
of n, around n w («) will have a noticeable contribution to the 
MDM. We also note that the average magnification from re- 
gions with n EIPs, (fi)„, is approximately linear in n (see Ta- 
ble 3), so that we expect A///(/x) <~ («)~ 1,/2 <C 1. Therefore, 
the width of the MDM is expected to be largest for (n) ~ 1. 
This seems to be in rough agreement with the results of nu- 
merical simulations (Wambsganss 1992; Irwin & Lewis 1995; 
Schechter & Wambsganss 2002). 
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5. APPENDIX A 

The inner integral in equation (6) may be solved analytically, 
reducing it to a single integral 



(p) 2tt 7_!_ 7 1 



A + BJ (nl+Sj) 



where (p) = |(1-k») 2 -7 2 | 1 [see equation (3)] and 
A = v /i + k2- 7 (2Si + 7) , B= Vl-^i + 7) 2 



(15) 
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This integral may be easily evaluated numerically. Alterna- 
tively, changing variables in equation (6) to r and <j> where 

x = 7 + Si, y = 52, r 2 =x 2 + y 2 and tan((fi) = y/x, gives 



1 r 2n 



dr 



r(\-r 2 ) 



(N) _ K* 

(fj) 2n J " J Q " r ( ttl + -f 2 + r 2 - 2 cos ffi/ 2 

2 K * (\ Kl-r 2 )EllipticE(f, 7 ^) 
= / dr — — =f- , 

* Jo [Kl + (r+j) 2 ]^Kl + (r-j) 2 



(17) 



where 



EllipticE(</> 



r<t> 

,x)= / d9(l-xsm 2 9) 1/2 . (18) 
Jo 



For the case of 7 = «*, corresponding to a singular isothermal 
sphere, we obtain an analytic result: 

= — !— (-6 7 2 EllipticE rarcsin(Q,C- 2 l 



-3(1- 7)7 2 EllipticF [arcsin(C) , CT 2 ] 
+ C[ 7 (27 2 -l)EllipticK(C 2 ) 

+ D EllipticPi(£ , C 2 ) + D*EllipticPi(E * , C 2 )] |l 9) 



where 
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1 +27+27 2 



D _ /+(l-07+(2-/)7 2 -(l + 307 3 



(20) 



l + (l + /)7 

D* (E*) is the complex conjugate of D (E), and 

r<t> 

EllipticF(0,x) = / d9(l-xsm 2 6y l/2 , 
Jo 

EllipticK(x) = EllipticF(7r/2,x) , 

j-TT/2 

EllipticPi(x,y) = / rf^l-xsin^-'a-y sin 2 6>r 1/2 (21) 
Jo 
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Table 1 
The ratio (N)/{h). 
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K*=0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1 


1.2 


1.5 





1 


0.8190 


0.6721 


0.5536 


0.4584 


0.3820 


0.3206 


0.2711 


0.2310 


0.1983 


0.1716 


0.1311 


0.0917 


0.1 


0.99 


0.8105 


0.6650 


0.5478 


0.4537 


0.3782 


0.3176 


0.2687 


0.2292 


0.1969 


0.1704 


0.1303 


0.0913 


0.2 


0.96 


0.7850 


0.6438 


0.5305 


0.4398 


0.3672 


0.3089 


0.2618 


0.2237 


0.1926 


0.1670 


0.1281 


0.0901 


0.3 


0.91 


0.7428 


0.6089 


0.5022 


0.4171 


0.3492 


0.2947 


0.2507 


0.2149 


0.1856 


0.1615 


0.1246 


0.0882 


0.4 


0.84 


0.6841 


0.5608 


0.4635 


0.3864 


0.3250 


0.2756 


0.2357 


0.2031 


0.1763 


0.1540 


0.1198 


0.0856 


0.5 


0.75 


0.6095 


0.5005 


0.4156 


0.3487 


0.2955 


0.2526 


0.2176 


0.1889 


0.1650 


0.1450 


0.1140 


0.0824 


0.6 


0.64 


0.5200 


0.4296 


0.3602 


0.3058 


0.2621 


0.2266 


0.1973 


0.1728 


0.1522 


0.1348 


0.1073 


0.0787 


0.7 


0.51 


0.4173 


0.3505 


0.2999 


0.2596 


0.2266 


0.1989 


0.1756 


0.1557 


0.1386 


0.1239 


0.1001 


0.0746 


0.8 


0.36 


0.3048 


0.2678 


0.2383 


0.2131 


0.1908 


0.1711 


0.1536 


0.1382 


0.1245 


0.1125 


0.0925 


0.0702 


0.9 


0.19 


0.1907 


0.1886 


0.1805 


0.1692 


0.1568 


0.1444 


0.1323 


0.1211 


0.1107 


0.1012 


0.0848 


0.0657 


1 





0.0974 


0.1236 


0.1314 


0.1310 


0.1265 


0.1200 


0.1126 


0.1049 


0.0974 


0.0902 


0.0771 


0.0611 


1.2 





0.0303 


0.0534 


0.0682 


0.0764 


0.0801 


0.0808 


0.0795 


0.0770 


0.0739 


0.0703 


0.0628 


0.0521 


1.5 





0.0111 


0.0212 


0.0298 


0.0366 


0.0415 


0.0448 


0.0468 


0.0477 


0.0478 


0.0472 


0.0449 


0.0400 





NOTE. — The ratio of the mean number of positive parity micro-images (N) and the mean total magnification (/x) for different 
values of the normalized surface mass density in stars, k*, and the external shear, 7. 7 and k*, as determined with equations (6), 
(15), (17), (19). 



Table 2 

Numerical Simulations Versus Theory 



K* =7 


( n )sim 


(«)th 


(M)sim 


G")th 


0.333 


0.367 


0.3806 


2.902 


2.994 


0.400 


0.901 


0.9320 


4.992 


5.000 


0.666 


0.531 


0.5759 


2.907 


3.012 



Note. — The first column labels the simulation, 
through the values of the normalized surface mass den- 
sity in stars, k», and the external shear, 7. The next 
two columns provide the average number of extra im- 
age pairs calculated from the simulation (according to 
equation 10), (n) s i m , and the theoretical prediction for 
this quantity, (n)th, calculated using equations (2) and 
(6). The last two columns are the average magnifica- 
tion of the simulation (fi) s i m , and the theoretical aver- 
age magnification (^)th (from equation 3). 
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Table 3 

Results of Numerical Simulations 



k* =7 


n 


fn 




Pn 









0.690800 


0.690800 


0.693511 


1.780 




1 


0.251230 


0.942030 


0.252216 


4.903 


0.333 


2 


0.049298 


0.991328 


0.049492 


6.999 




3 


0.003878 


0.995206 


0.003893 


11.62 




4 


0.000777 


0.995983 


0.000780 


15.07 




5 


0.000108 


0.996091 


0.000108 


17.14 







0.362692 


0.362692 


0.367695 


2.387 




1 


0.412717 


0.775409 


0.418409 


5.352 




2 


0.164497 


0.939906 


0.166766 


7.726 


0.400 


3 


0.039863 


0.979769 


0.040413 


10.51 




4 


0.005832 


0.985601 


0.005912 


13.16 




5 


0.000749 


0.986350 


0.000759 


16.65 




6 


0.000045 


0.986395 


0.000046 


26.33 







0.584765 


0.584765 


0.589852 


1.239 




1 


0.300253 


0.885018 


0.302865 


4.497 


0.666 


2 


0.093123 


0.978141 


0.093933 


6.807 




3 


0.012753 


0.990894 


0.012864 


9.935 




4 


0.000482 


0.991376 


0.000486 


14.90 



Note. — The first column labels the simulation, through the 
values of = 7. The remaining five columns are the number of 
extra image pairs n, the corresponding fraction of the source plane 
/„, its cumulative value up to n, its normalized value p„ calculated 
according to equation (9), and the mean magnification of regions 
with n extra micro-image pairs, ( / u)„ (that are marked by inverted 
triangles in Figure 4). 
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10 



k* = y = 0.333 





1 2 3 4 5 



K* = y = 0.400 




2 3 4 5 6 




1 2 3 4 5 6 



k* = y = 0.666 




n 




n 



FIG. 3 . — The bars show the probability p„ of having n extra image pairs, calculated from our simulations (k* = 7 = 0.333, 0.400, 0.666), using equation (9). The 
lines correspond to a Poisson distribution [equation (13)] where the solid line is for the analytical (theoretical) mean value (n), and the dashed line is for the mean 
value from the simulation («) s j m , calculated according to equation (10). The theoretical values of p„ for integer values of n are shown by the asterisk and circle 
symbols. The panels on the left (right) column are linear (logarithmic) representations of p„. 
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FIG. 4. — A log-log plot of the magnification distribution for the macro-image fip(p) (uppermost line), for our simulations, and its decomposition, /j,p n (/j,), to the 
contributions from regions with different numbers of extra image pairs, n. The alternating thin and thick solid lines, from left to right, are for n = through 6. The 
inverted triangles represent the average magnification from regions with n extra image pairs, (/i)„. 



